Numerical simulations of string networks in the Abelian-Higgs model 
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We present the results of a field theory simulation of networks of strings in the Abelian Higgs model. 
Starting from a random initial configuration we show that the resulting vortex tangle approaches a 
self-similar regime in which the length density of lines of zeros of <j> reduces as t~ 2 . We demonstrate 
that the network loses energy directly into scalar and gauge radiation. These results support a 
recent claim that particle production, and not gravitational radiation, is the dominant energy loss 
mechanism for cosmic strings. This means that cosmic strings in Grand Unified Theories are severely 
constrained by high energy cosmic ray fluxes: either they are ruled out, or an implausibly small 
fraction of their energy ends up in quarks and leptons. 
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In recent years, a sustained effort has gone into under- 
standing the formation and evolution of networks of cos- 
mic string, principally to provide a mechanism for seeding 
gravitational collapse [QJ|] . 

The picture that emerges from lattice simulations of 
string formation is a network consisting of a small number 
of horizon crossing self avoiding random walks, together 
with a scale invariant distribution of loops || . The sub- 
sequent evolution is driven by a tension in the strings 
causing them to straighten out. When two lengths of 
string pass through each other, they may intercommute, 
that is exchange partners. This allows the production of 
loops which can decay through gravitational radiation or 
particle production, depending on their size. 

A number of numerical studies of network evolution 
have been carried out using the Nambu-Goto approxi- 
mation of the string as a 1-dimensional object with a 
tension \i. In particular we will refer later to work by 
Allen and Shellard (AS) and Bennett and Bouchet (BB), 
both reported in 1. A consensus emerged from these 
studies that the network will relax into a scaling regime, 
where the large scale features of the network (the inter- 
string distance £ p and the step length £ p ) grow with the 
horizon in proportion to t. On scales less than £ p , there 
was evidence for a fractal substructure covering a range 
from the resolution scale to £ p , which seems to be the re- 
sult of kinks left behind by the production of loops. The 
loop production function itself did not scale: the distri- 
bution of loop sizes reported by BB is peaked around 
a cut-off introduced by hand into the simulation to en- 
sure that even the smallest loops are approximated by a 
reasonable number of points. 

It has been thought that the build-up of small scale 
structure will only be stopped by back-reaction from the 
string's own gravitational field when the intermediate 
fractal extends to scales of order G/x£ p , where G is New- 
ton's constant. Then the small scale structure and loop 



production will scale with the horizon. Although much 
smaller than £ p , loops produced at this scale are still 
vastly bigger than the string core and will decay through 
gravitational radiation. 

This picture has some support from a detailed analyt- 
ical study by Austin, Copeland and Kibble || (ACK) 
which deals in part with length scales £ p , £ p and a scale 
C, which can be interpreted as an angle-weighted average 
distance between kinks on the string. However, the ACK 
analysis relies on a number of unknown parameters, and 
for some parameter ranges small scale structure is absent. 

In ||, Vincent, Hindmarsh and Sakellariadou (VHS) 
suggested a different picture based on results from 
Minkowski space Nambu-Goto simulations using the 
Smith- Vilenkin algorithm |7j . This algorithm is exact for 
string points defined on a lattice, allowing easy detection 
of intercommutation events. When loop production is 
unrestricted (up to the lattice spacing) they found that 
small scale structure disappeared, although loop produc- 
tion continued to occur at the lattice spacing. They could 
only recover the small scale structure seen in (J] by arti- 
ficially restricting loop production with a minimum loop 
size greater than the lattice spacing. The suggestion is 
that the small scale structure seen in other simulations is 
an artifact of a minimum loop size and that loop produc- 
tion will occur at the smallest physical scale: the string 
width. 

This has interesting consequences for energy loss from 
string networks. Loops formed at this size will decay into 
particles rather than gravitational radiation and could 
give a detectable flux |0] . If particle production provides 
the dominant energy loss channel, then GUT theories 
with strings are heavily constrained. 

In this letter, we present the results of a series of field 
theory simulations of networks of Abelian Higgs vortices. 
They support the results of Vincent et al . If loops form, 
they form at the string width scale and promptly col- 
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lapse. However, most of the string energy goes directly 
into oscillations of the held - radiation. We measure lit- 
tle small scale structure in the network. Furthermore, we 
find that the network scales with a scaling density consis- 
tent with the Smith- Vilenkin simulations, implying that 
the latter simulations do not unduly exaggerate energy 
loss as a lattice effect, a possibility suggested in ||. 

The simplest gauge strings are contained in the 
Abelian-Higgs model, which has a Lagrangian 

C = (2V) f W) ' ~ ^(H 2 - * 2 ) 2 , (1) 

where is a complex scalar field, gauged by a U(l) vec- 
tor potential with a covariant derivative = — ieA^. 

is the field strength tensor <9mA„ — d v A^. To 
model the system described by Eq. ffl we use techniques 
from hamiltonian lattice gauge theory Q (and references 
therein) . 

An attractive feature of this formalism is that the 
Hamiltonian respects the discrete version of the U(l) 
gauge transformations and consequently if Gauss's law 
is true initially, then it is true for all later times. Signifi- 
cant violations of Gauss's law are obtained if instead one 
uses a numerical scheme based on finite-differencing the 
Euler-Lagrange equations for <f> and Aj. 

The equations of motion derived from the Hamiltonian 
are evolved by discretising time with At = ra (r -C 1) 
and using the leapfrog method of updating field values 
on even time steps and the conjugate momentum on odd 
time steps. In keeping with the systems we are trying 
to model, we create initial conditions by allowing an en- 
ergetic configuration of fields to dissipate energy until cf> 
is close to the vacuum everywhere except near the string 
cores. We are not attempting to model the formation 
process, rather we wish to create a reasonable random 
network of flux vortices for subsequent evolution. 

The simulation proceeds as follows. We generate a 
Gaussian random fc-space configuration of each compo- 
nent of <j) with a /c-dependent Gaussian probability dis- 
tribution with variance Ok = \Jmbj (k 2 + m 2 ). Sam- 
ples drawn from this distribution are then Fourier trans- 
formed to x-space to generate a starting configuration, 
which will be uncorrelated on scales greater than m _1 . 
By varying the two paramenters m and b it is possible 
to control both the average amplitude and the correla- 
tion length of the initial configuration (this is in turn 
related to the initial defect density after the dissipation 
period). This then forms the initial conditions to be dis- 
sipated along with <p = & 1 = E l = 0. Our dissipation 
scheme is to add the gauge invariant terms rjir x and r\E x 
to the equations of motion for <f> and 9 l . This ensures 
that Gauss's law is preserved by the dissipation, but 
changes the charge-current continuity equation so that 
p + V • j = —rjp. We have found that the most efficient 
dissipation is when the charge grows from to e during 
the dissipation process. This allows the early formation 
of essentially global vortices, into which the magnetic flux 
relaxes as e is increased. 



Once the network has formed with (\4>\) « 1, we evolve 
the network with the dissipation turned off for half the 
box light-crossing time. During the evolution energy is 
conserved to within sw 0.3%. 

A reasonable resolution of the string core is obviously 
an important requirement of the simulation. Tradition- 
ally the core is defined by the inverse masses of the field, 
so for e = 1, A = 2 string, d c — 1. However, the fields de- 
part appreciably from the vacuum over a larger distance, 
about 4 units. 

We have performed a series of simulations on lattices 
of 224 3 , 304 3 and 336 3 points, lattice spacings of between 
0.25 and 0.75, and with r = 0.15. We vary the parameter 
m in the variance of the field distribution to control initial 
defect density. In this letter we keep A/2e 2 = 1. 

Due to the size of the simulations we are unable to 
achieve the same level of statistical significance obtained 
using Nambu-Goto codes, although there are some re- 
vealing qualitative results. 
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FIG. 1. Plots of £ p for a series of 336 simulations with 
different lattice spacings. From top to botttom a —0.75, 0.65, 
0.75, 0.4, 0.5, 0.45 and 0.25. In the initial gaussian distri- 
bution the parameter m is varied to give the different initial 
correlation lengths. £ is given in units of the inverse scalar 
mass m^ 1 . 

The gauge-invariant zero of the <fi field provides a sim- 
ple way of analysing the network of vortices. A string 
passes through a lattice cell face if the winding around 
the four corners is non-zero. By starting in any lattice cell 
with a string in it, the string can be followed around the 
lattice until it returns to the starting point. The distri- 
bution of box-crossing string and loops can then be anal- 
ysed in much the same way as for Nambu-Goto strings 
PjQ< In particular, we are interested in the behaviour 
of the length scale £ p , roughly the inter-string distance, 
defined as £ p = V/i| , where L^ p is the total physical 
length of string above length £ p (not the more usual in- 
variant length). In practise, we calculate the length of 
string by smoothing over the string core width d c . Scal- 
ing occurs if £ p grows with t, although the actual value of 
x p in the scaling density p^ p = /i/£ 2 = p/x^t 2 , depends 
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on the efficiency of the energy loss mechanism. Fig. |l| 
shows the function £ p for a sample set of runs with dif- 
ferent simulation parameters. Although with small dy- 
namic ranges, one can never be sure that £ p is approach- 
ing scaling or just on a slow transient, the scaling values 
for x all appear to be in the range 0.27 — 0.34, which is 
in agreement with Smith- Vilenkin simulations with max- 
imum loop production where x p = 0.27 ± 0.05. This 
value is obtained by converting the invariant length 
used to define £ in || to physical length L^ p , if X'(cr) 
is the string tangent vector, then the physical length is 
just = J do"|X'|. In the latter half of the simulation 
runs, these plots are extremely linear: the exponent of t 
in £ oc t v isp = 1.002 ± 0.03. 

In H , VHS argued that the string network is be dom- 
inated by processes on the smallest physical scale, ei- 
ther radiation or the production of loops at the string 
width. In all our runs performed so far, we find that 
this is true. When analysing the distribution of string 
formed by zeros of <j>, we consider three classes: long 
string [L > £ p ), intermediate loops {%d c < L < £ p ) and 
core loops (L < Trd c ). Fig. shows a typical run for with 
a proportionate break down into types of loop over time. 
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FIG. 2. Proportion of string length intermediate loops 
(grey) and core loops (black) . The long string (white) contri- 
bution makes up 100%. The lattice spacing is a = 0.35. 

Over the first part of the simulation we see an overall 
decrease in the percentage of string in loops of both types. 
Over the latter part of the simulation we see three sharp 
increases in % of intermediate loops at t = 37, t = 61 
and t — 81. In all three cases, this is caused by a single 
large loop (with L > £ p ) decaying to below the £ p thresh- 
old. Although these events constitutes a loss to the long 
string network, the decay channel is radiation and not an 
intercommutation event producing an intermediate loop. 
Note that aside from these loops, in the latter half of 
the simulation, all the string length is in long string (to- 
gether with a small population of core loops), and the 
network continues to scale, despite some minor fluctua- 
tions in £ p corresponding to the shrinking large loops (see 
curve marked 'x' in Fig. |l|). We infer that all energy loss 



occurs ultimately through radiation and that it is very 
efficient at scaling the network. 

To back up our claim that radiative processes are key 
to network evolution, we examined sinusoidal standing 
waves with wavelength A and initial amplitude A where 
A = A/ 2. The necessary scaling behaviour for the en- 
ergy density in a network can be seen in the radiation 
from this system (incidentally one without intercommu- 
tation). Consider that a string network can be thought 
of as a length £ in a volume £ 3 , giving a density p = /i/£ 2 
and p = — 2^£/£ 3 . For a scaling network, p£ 3 = con- 
stant. We studied a series of standing waves with in- 
creasing wavlength A and measured the energy density 
p{t) in a sub-volume V\ ~ A 3 around the string. We 
found that p decays fairly linearly as the radiation from 
the string begins to leave V\, and then tails off as A 
decreases significantly. We measure p during the linear 
region when A s» A/2. We find that p\ 3 is constant to 
within 10 percent over a range of standing waves with 
A from 20 to 180, two orders of magnitude bigger than 
the string width. We also note that the pattern of oscil- 
lations around the standing wave remains substantially 
unchanged by an increase in the lattice spacing: we be- 
lieve that any differences can be accounted for by the 
resulting decrease in resolution. The fact that the scal- 
ing density in the network simulation does not depend 
on the lattice spacing is also evidence that the radiative 
process is not a lattice artifact. 
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FIG. 3. The fractal dimension of the long string at two 
times late in the evolution for B = 336 3 and a = 0.75 (solid), 
a — 0.25 (dashed). The curves do not coincide, which may be 
due to difference in dynamic range in the two simulations. 

It is claimed in ^] that if loop production is unre- 
stricted, small scale structure will disappear. Small scale 
structure has been studied in two ways. Firstly, there is 
a length scale £ mentioned in the introduction. Secondly, 
the relation between end-to-end distance R and distance 
along the string D, (D/R) oc (Z?/£ p )™, defines a scale- 
dependent fractal dimension n. If strings are smooth on 
small scales and brownian on large scales, n should in- 
terpolate from ~ to ~ 0.5. BB and AS observed that 
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over a significant fraction of the range from the resolution 
scale to £, n is constant. VHS showed from the behaviour 
of £ that the region of constant n did reflect the existence 
of small scale structure: when loop production is unre- 
stricted the region of constant n is absent and £ ~ £; 
when loop production is restricted a region of constant n 
emerges and £ -C £. 

Fig. H shows that for the flux vortices n varies smoothly 
from ~ 0.01 to ~ 0.5. The absence of the transient is in 
good qualitative agreement with the Nambu-Goto sim- 
ulations with £ ~ £, and we can argue from this that 
gravitational radiation (gr) is irrelevant to our claimed 
particle production (pp). Energy loss processes can be 
described by an efficiency parameter c in the rate equa- 
tion p = —cp/£. For a scaling network with £ = e£, 
c gr = TG^/e where T - 10 1 - 10 2 . If there is no 
small scale structure, e ~ 1 and, for GUT strings, 
c gr = 10~ 4 — 10~ 5 . This is to be compared with the 
overall energy loss through radiation in our simulations 
which gives c pp = 0.6. It is only when there is consider- 
able small scale structure (e ~ 10 -4 ) that gravitational 
radiation becomes significant. 

We believe that our results support the claim of VHS 
that the dominant energy loss mechanism of a cosmic 
string network is particle production. Our field theory 
simulations show that the string network approaches self- 
similarity, as £ p oc t, with a constant of proportionality x p 
consistent with Nambu-Goto codes which simulate string 
objects directly. We provide evidence that the long string 
network loses energy directly into oscillations of the field, 
which is the classical counterpart of particle production, 
or into loops of order the string width which quickly 
decay into particles. Previous calculations on particle 
production underestimated its importance by using per- 
turbative calculations with light particles p0[ . We see 
non-perturbative radiation at the string energy scale. 

This represents a radical divergence from the tradi- 
tional cosmic string scenario, which held that the loops 
are produced at scales much greater than the string 
width and subsequently decay by gravitational radia- 
tion. It implies that strong constraints come from the 
flux of ultra high e nerg y (UHE) cosmic rays. Bhat- 
tacharjee and Rana estimated that no more than 
about 10~ 3 of the energy of a GUT-scale string net- 
work can be injected as grand unified X particles (as- 
sumed to decay into quarks and leptons). More recent 
calculations [O] limit the current energy injection rate 
Qo of particles of mass mx to approximately Qq < 
2 x 10- 34 (m x /10 15 GeV) - 44 ergss- 1 cm- 3 (this formula 
is based on an eyeball fit to their most conservative bound 
in Fig. 2 Jn|). Sigl etal |ll| obtain somewhat lower val- 
ues for Qq: we have been cautious and adopted the higher 
value as our bound. 

The energy injection rate from a scaling network is of 
order [i/t 3 ~ 10" 32 ergs s^cm" 3 (/x/10 32 GeV 2 ). Tak- 

-29 f -1.3 Ge y2 Qr 



ergy appearing as quarks and leptons. Thus we conclude 
that GUT scale strings with Gp ~ 10~ 6 are ruled out: 
we regard a fraction /x~5x 10~ 3 as implausible. 
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10m x , we derive \i < 10 
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G\i < 10 9 fv 1 ' 3 , where fx is the fraction of the en- 
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